Jahn- Teller induced Berry phase in spin-orbit coupled Bose-Einstein condensates 
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We demonstrate that Berry phases may greatly affect the dynamics of spin-orbit coupled Bose- 
Einstein condensates. The effective model Hamiltonian under consideration is shown to be equivalent 
to the E X e Jahn- Teller model first introduced in molecular physics. The corresponding conical 
intersection is identified and the Berry phase acquired for a wave packet encircling the intersection 
studied. It is found that this phase manifests itself in the density profile of the condensate, making 
it a directly measurable quantity via time-of-flight detection. Moreover, the non-Abelian gauge 
structure of the system is addressed and we verify how it affects the dynamics. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensates flj 
opened a new door for the study of quantum phenom- 
ena. Since then the field has seen a tremendous progress 
in terms of trapping, cooling and manipulating atoms, 
covering cold atoms in optical lattices [2], BCS-BEC 
crossover |3] and multi-component Bose-Einstein con- 
densates [4]. Moreover, it has been demonstrated both 
theoretically 5, 6] and experimentally 0] that ultracold 
multi-level atoms moving in spatially dependent laser 
fields give rise to effective gauge potentials. Consider- 
ing four-level atoms in a tripod setup, non-Abelian gauge 
fields are achieved [1] . The idea of utilizing a four-level 
atom with three degenerate ground states coupled to one 
excited state to generate non-Abelian gauge potentials 
dates back to the work of Unanyan and co-workers @, 
where however the effective magnetic fields arise from 
time-dependence of the Hamiltonian rather than spa- 
tial dependence. The particular system configuration of 
Ref. [8,] has turned out to be extremely rich, and a se- 
ries of papers investigating various aspects of the model 
have been published recently. Among these are, spin- 
Hall effects [13], the Aharonov-Bohm effect rela- 
tivistic characteristics [l^, [I^! spin-echo phenomena of 
trapped fermions spin dynamics [l5| . expansion to 
time-dependent laser fields [1^1, novel phases of the con- 
densate [13], the structure of the energy spectrum p^ . 
and effective magnetic monopoles [T9| . 

For certain laser arrangements, the energy dispersions 
possess a point- degeneracy also termed conical intersec- 
tion (CI) po| . It is well known [U [S^l that apart from 
the dynamical phase, encircling a CI renders a Berry 
phase of the state vector Differing from typical 

situations in molecular and chemical physics, here the 
CI occurs in momentum and not in position space. As a 
consequence, the Berry phase is more directly manifested 
in the momentum wave function. Similar situations ap- 
pear in graphene [131 and cold atoms in elaborate optical 
lattice configurations [25| . 



In the current work, we consider a harmonically 
trapped spinor Bose-Einstein condensate in the presence 
of spatially dependent laser fields such that a CI is re- 
covered in momentum space. By taking the harmonic 
confinement of the condensate into account, the result- 
ing Hamiltonian is equivalent to the E x e Jahn- Teller 
one [1^ , but with position and momentum interchanged. 
Moreover, the non-linear terms arising from atom-atom 
scattering are included. Such terms are absent in for ex- 
ample CTmodels in molecular physics [l^, [13, [11] . Uti- 
lizing numerical wave packet simulations, we thoroughly 
discuss the Berry phase and how it affects the system 
dynamics, both on a short and long time scale. It is 
found that over longer time periods, the population of 
the two phonon modes is swapped, and in particular, the 
characteristics of this exchange mechanism depend on 
the Berry phase. Various studies have been concerned 
about Berry phases in Bose-Einstein condensates influ- 
enced by slowly varying external fields [1^, [13] . Con- 
trary to these references, in the present model the Berry 
phase shows up in the internal spinor dynamics as an ef- 
fective spatially dependent gauge field, in analogy with 
the Born-Oppenheimer scenario in molecular theory [2]| . 
We note that Berry phases were as well considered in 
Refs. [lll[3, but without taking into account for a trap- 
ping potential nor atom-atom collisions. A harmonic 
trapping potential causes the momentum wave packet 
to spread, which has important consequences for its dy- 
namics. In addition to investigating the consequences of 
a non-zero Berry phase, we also study the non-Abelian 
structure of the system. We show that one round-trip of 
the CI clockwise or anti-clockwise brings about different 
final states. 

The paper is outlined as follows. In the next sec- 
tion, the model system is introduced and the effective 
Gross-Pitaevskii equation given. The equivalence with 
the E X £ Jahn- Teller Hamiltonian is demonstrated and 
we discuss the form of the momentum potential energy 
surfaces. Section Hill is devoted to numerical simulations. 
First, in Sec. IIII Al the short time evolution is consid- 
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ered, the non-cyclic Berry phase [3l| for the spinor state 
is computed, the impact of the Berry phase on the wave 
packet dynamics is clarified, and we show how a time- 
of-flight measurement of the condensate would reveal the 
presence of a Berry phase. The following Sec. IIII 51 deals 
with the long time scale properties. A collapse-revival as 
well as a swapping phenomena is found. The underlying 
non-Abelian structure is envisaged in Sec. IIII CI Last we 
finish with concluding remarks in Sec. IIVI 



II. MODEL SYSTEM 

We begin by deriving the effective single atom Hamil- 
tonian, and then turn to the many-body case by includ- 
ing atom-atom scattering to obtain a mean field Gross- 
Pitaevskii equation. 

The setup is detailed in Fig. [1] (a). Three degenerate 
meta stable Zeeman ground states are dipole coupled to 
an excited state via three respective external lasers. Us- 
ing the notations of the figure and imposing the rotating 
wave approximation, the interaction Hamiltonian reads 



Hi = hA\o){o\ + h{ni\o){i\ 

+r!3|0)(3|-f h.c). 



r!2|0)(2| 



(1) 



which defines the full Hamiltonian H — p'^ /{2m) -f i//, 
p and TO being the momentum and mass, respectively, 
of the atom. The coupling strengths, being propor- 
tional to the laser amplitudes, are parameterized as 
fii = sin 6* cos (/9e''^i , Vl^ — sin sin (/?e''^^ and = 
n cos 6*6*^3 with n = ^|f]i|2 + |f]2p + |f^3p- Here, 6*, ^ 
and Si are allowed to be spatially dependent. The in- 
teraction part Hi of the full Hamiltonian may be readily 
diagonalized, rendering the eigenstates \wq), I^^s): I""!); 
and |m2), where in particular the two states 
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are degenerate with zero eigenvalue, and here 5"^ — 
Si — Sj . These two states are termed dark states as they 
do not couple to the excited state |0) via Hi. Due to 
the shifted energy of the remaining two states \wo) and 
juia) compared to and |m2), they can be adiabatically 
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FIG. 1: Atomic (a) and laser (b) configurations. 



eliminated from the Hamiltonian provided that A ;:!> f2, 
ending up with an effective two- level problem for \ui) and 

l"2). 

Before writing down the resulting Hamiltonian we 
specify the laser configurations Jjj, see Fig.[T](b). First 
we assume Si = iS'2, 5*31 = rnvgy/h and ip = mv^px/h 
and 9 € [0, tt] some constant. Two lasers, which govern 
the coupling of |1) and |2) to |0), are propagating along 

ki and k2 in the iy-plane, where the angle of intersec- 
tion between them is C. In the i-direction, these two 
lasers are standing waves, however mutually shifted in 
phase by 7r/2, and propagating waves in the y-direction. 
A third laser is a propagating wave along fca in the yz- 
plane, whose direction is determined by the angle ^ rel- 
ative the y-axis. In terms of the angles C a-nd ^ and the 
wave numbers ki = I A:,- 1 we have 



^i = 52 = fciycos(C/2), 
<5'3 = ^32; cos ^, 
if = 2kix sin(C/2). 



(3) 



After the adiabatic elimination, the effective single par- 
ticle two-level Hamiltonian becomes 



HcS = T: 1- So^y 

2m 



voCTxPx - viayPy, 



(4) 



where Sq is an effective Zeeman splitting, vq — cos 0, 
WsSin^^/2 and (t^, with i = x,y,z are the Pauli 

,y _ -i\ui){u2\ ' 



Vi 



x,y,z 

operators a^; = |uD('"2| + 1^2) 
«|m2) and (7z — \ui){ui\ — \u2){u2\- The p-dependence 
of the last two terms originates from non-adiabatic cou- 
plings [13], and we note that these do not include p^ 
and we may therefore factor out the ^-dependence. The 
spin-orbit coupling is on a Rashba-ioim [33], which has 
been frequently studied in the context of semiconductor 
quantum-dots [s^-lsH- 

We will further add a harmonic trapping potential (as- 
sociated with angular frequency lo) for the atoms. This 
was in fact also done in Ref. jS^]- However, for the con- 
densate considered in this work, such a trapping poten- 
tial is almost exclusively present in experiments, while 
this is not the case for the electrons in semiconducting 
quantum-dots of Ref. [35| . By properly Stark shifting 



the internal atomic states, individual trapping poten- 
tials can be achieved. In particular, for certain constant 
Stark shifts, the Zeeman splitting term in the Hamilto- 
nian [TT, T^l can be made to vanish, which we will assume 
in the following. It should be mentioned that in order to 
justify adiabatic elimination of the two dark states, once 
the harmonic trapping potential has been included, one 
should assume not only A 3> il, but also ^ uj in order 
to ensure the validity of the Born-Oppenheimer separa- 
tion. 

Imposing s-wave scattering between the atoms in terms 
of the regular non-linear term gives a Gross-Pitaevskii 
equation. Defining the spinor wave function 



(i,yi(uii*(i)y 

(i,y|(u2|*(t)). 



(5) 



where tpiljcjij^t) — {x , y\Tpi{t)) , i ~ 1, 2 and the nor- 
malization reads J dxdy\'^{x,y, = 1, the Gross- 
Pitaevskii equation takes the form 



ihM{x,y,t)-- 



2m \ dx 



\di^ ^ W ) 



dx 



ihvicr 



Vdy 



+^Hx,y,t)g^ix,y,t^^{x,y,t), 

(6) 

where cr^ are the standard Pauli matrices acting on the 
spinor space and 



511 512 

512 522 



(7) 



is the matrix valued scattering amplitude. In general, 
the scattering amplitudes gu, ^22 and 512 can be differ- 
ent. We will in the following assume 511 = 322 = 5 
and gi2 = 0, a condition that is supposed to be approx- 
imately achievable for a ^^Rb condensate [sj]. We have 
numerically verified, using a non-zero gi2 does not give 
rise to any qualitative changes of our results. In terms 
of physical quantities, g = Airfi^Na/mV, where N is the 
number of condensate atoms, V the effective volume and 
a the s-wave scattering length. Attractive atom-atom in- 
teraction is obtained for ^ < and repulsive for g > 0. 
Noteworthy is that the parameter g may be tuned over a 
wide range of values using Feshbach resonance techniques 

The single particle Hamiltonian Eq. Q, once the har- 
monic trapping potential has been included, is equiva- 
lent to the E X e Jahn- Teller one, which was first intro- 
duced in molecular physics [26j . However, the canoni- 
cal momentum and position variables are interchanged 
in the original E x e Hamiltonian. Mathematically this 
is just a question of representations, but it has physical 
consequences, for example in state preparation and mea- 
surement. Another outcome of this is that the effect of 
the Berry phase normally apparent in the physical state 
'i'{x,y,t) will instead manifest itself in the spinor wave 



V(p) 




FIG. 2: Effective APS corresponding to the Gross-Pitaevskii 
equation (O. The conical intersection is located at the origin, 
Px = Py = 0, where the two APS become degenerate. Here 
we have vo = vi. 



function in the momentum representation, $ 

being related to ^ via the Fourier transform 

HpB,Py,t)= [ dxdye-'^''P'+yPy^/''^{£,y,t). (8) 



The phenomena we are studying are therefore more con- 
veniently extracted by analyzing the problem in a mo- 
mentum representation, and in particular in terms of an 
effective momentum-dependent potential. 

For this purpose we introduce the unitary adiabatic 
transformation [/ad — Us.d{ps:,Py) that diagonalizes the 
spin-orbit term, namely 



ad 



VoPi - iVlPy 

vqPx + iviPy 



U. 



(9) 



= \J (voPif + (viPyY 



EE Act,. 



This results in the effective adiabatic potential surfaces 
(APS) in momentum space, defined as 



ad 



1 (^2 



± A 



(10) 



which conically intersect at p = 0. Encircling the conical 
intersection (CI) gives rise to a Berry phase [1^, which, 
indeed, is related to the gauge properties of the problem 
t21i, |28(]. Figure displays the APS in the polar symmet- 
ric case in which vq = vi. For vq ^ vi, on the other hand, 
the polar symmetry is broken and the lowest APS pos- 
sesses two local minima instead of the symmetric Mex- 
ican hat sombrero shape. In what follows, we will only 
consider the situation with vq = vi = v, a condition that 
can be met by adjusting the laser parameters. Thus, in 
our simulations, the lower APS in momentum space has 
the familiar sombrero structure. 
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III. DYNAMICS 

We numerically solve the Gross-Pitaevskii equation ([6]) 
using the split-operator method [s^l- First, we intro- 
duce dimensionless variables scaled by the characteristic 
length and energy given by / = yjh/mui and Ec — huo, 
respectively. The Gross-Pitaevskii equation then attains 
the form 



to momentum space via the fast Fourier algorithm and 
the second exponent multiplies the momentum wave 
packet. Finally, the inverse Fourier transform gives the 
propagated state '^{x,y,T + St). The time step St is 
chosen such that contribution from the non-commuting 
part, due to factorizing the time-evolution operator, can 
be safely neglected. 



z|:*(x,y,r) = 



1 fsL 

2 \dx^ 



' dx 



d_ 

dy 



id_ _ d_ 

'' dx dy 





-h7(|V'i(a:, y, T)|2+|?/'2(a;, y, t)|2)^(x, ?/,r). 



where 



I ' 



T — tuj. 



(11) 



(12) 



a 



As initial spinor components, we take minimum uncer- 
tainty Gaussians 



V'i(a;,j/,0) = fli 



1 



1/2 



-i{Pxox+Pyoy) 



(x-xo)'' + (i/-i/o)^ 



xe 



(13) 



where |aip, normalized as |aip-|-|a2p = 1, determine the 
initial population in states \ui) (i = 1, 2), A/ is the initial 
wave packet width and xq, yo, Pxo and Pyo are initial 
averages of positions and momentum. We will choose 
A/ = 1, corresponding to the ground state width of the 
oscillator potential ^{x'^ + u'^)- 

An advantage of utilizing the split-operator method is 
that the time-dependent solution of the GP equation is 
automatically obtained in both momentum and position 
representation. As already pointed out, the effect of the 
Berry phase is most easily extracted from the momentum 
wave packets. More precisely, the split-operator method 
makes use of the fact that for sufficiently short time steps 
St, the time-evolution operator can be factorized into one 
momentum-dependent and one spatially dependent part 



\E'(a;, y, t+St) = cxp<— z 



xexp 



1/ 2 



-Py)+VPx(Tx+VPy(Ty 



-{x^ + y^) + :i-^{x,y,T)\^^ 



St 



(14) 

in the limit St 0. The first exponent multiplies 
\E'(x,?/,r), then the resulting wave packet is transformed 



A. Berry phase over short time periods 

Let us first examine how the Berry phase of the spinor 
part of the state behaves when the wave packet encircles 
the origin {px,Py) = (0,0). We mainly consider the case 
with no atom-atom scattering (g — 0), and only briefly 
discuss the g case. The spinor state is given by the 
reduced density operator 



p{t) = Pi{t)\ui)(ui\+P2{t)\u2){u2\ 

+C{t)\ui){u2\+C*{t)\u2){ui\ 



(15) 



where P,(t) = (V'»(t)|?A,(t)), i = 1,2, and C{t) = 
{ipi(T)\'ip2iT)) represent the populations and interfer- 
ence, respectively, of the two dark states ui,U2- By 
writing p(t) on Bloch form, we identify the time de- 
pendent Bloch vector r(T) = (m(t), u(r), w(r)) = 
(2ReC(r), -2ImC(T), Pi (r) - P2{t)). We choose an ini- 
tial wave packet characterized by the parameters ai — 
-02 = 1/V2, xo = 0, yo = -3, Pxo = v = 4, and Pyo = 0. 
For this choice, the wave packet will predominantly pop- 
ulate the lower APS and in particular sin ce P xo = w it is 
located at the minima of the sombrero [40,|4l[. Moreover, 
due to the non-zero yo < 0, as time evolves it traverses 
the sombrero minima clockwise. The motion of the Bloch 
vector representing the spinor state is shown in the upper 
panel of Fig. ([3]). Note that |r(0)| — 1 (pure state) since 
|f/;i(0)) = \'ip2{0)) = 1-0(0)), which implies that the initial 
spinor wave function takes the product form |4'(0)) = 
IV'(O)) (|ui) - \U2))/V2. The Jahn-Teller coupling cre- 
ates entanglement between the spinor and spatial degrees 
of freedom when the wave packet evolves, which explains 
why |r(r)| varies. To compute the evolving Berry phase 
7(r) for this effective non-unitary evolution of the spinor 
subsystem we write the reduced density operator on spec- 
tral form p(t) = Xk{T)\(pk{T)) {(pk(T)\ and use [3i| 



7(t) = arg 



J2V^kmkir){ipkmMr)) 



xe 



-Io{Vk{T')\Vk{T'))dr 



arg[(^i(0)|^i(r))] 

r{Mr')\viir'))dT', 
Jo 



(16) 



where the second equality follows from that p(0) is pure 
(only one Afe(O) is nonzero) and the assumption that 
Ai(r) is the largest eigenvalue over the relevant time in- 
terval. Thus, the Berry phase in this case becomes the 
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FIG. 3: Bloch vector and corresponding non-cyclic Berry 
phase for the spinor subsystem with no atom-atom scattering 
{g — 0). The initial state is characterized by the dimension- 
less parameters ai = — a2 = l/\/2, xo — 0, yo ~ —3, pxo = 4, 
and Pyo = 0. 



standard non-cyclic geometric phase [3l| for the eigen- 
state of p{t) with largest eigenvalue. Here, this phase 
becomes — ^ times the solid angle il,gc{T) enclosed by 
the path r' i— > n(T') — r(T')/|r(T')| and the shortest 
geodesies that connects its end-points n(0) and n(T). 
From the upper panel of Fig. Owe see that ^^gc(''') ~ 
for the first half of the path, and r2gc(r) « 27r for the 
second half. Thus, we expect 7(t) = — irigc(T) to make 
a TT-jump at half the period of the motion. This behavior 
of the non-cyclic Berry phase in the E ® e Jahn- Teller 
system has been examined in the literature 39'] and is 
verified in the lower panel of Fig. [3l Note the wiggles in 
the vicinity of the phase jump which are caused by small 
non-adiabatic effects in the wave packet motion. 

To extract measurable effects of the non-zero Berry 
phase we need the wave packet to self-interfere. To this 
end, we assume as above ai = —a2 = l/'s/Z, xq = yo = 
PyO = 0, Pxo = V, while now yo = 0. Thus, the mo- 
mentum wave packet has no initial velocity, so for zero 
atom-atom scattering corresponding to g = 0, the wave 
packet starts to spread around the minima of the som- 
brero. We emphasize that this spreading takes place in 
the momentum wave packet $, and is driven by the un- 



certainty in the spatial coordinates. Thus, such spreading 
is only possible since we consider a harmonically trapped 
condensate, contrary to Refs. [Ill where no trapping 
was taken into account. For positive g, the interaction 
is repulsive which tends to increase the spreading of the 
^ff wave packet while slowing down the <i> wave packet 
broadening. The opposite occurs for attractive interac- 
tion. Consequently, no overall spreading in $ takes place 
if g is large and negative. 

We assume a g such that for sufficiently long times, <& 
has spread noticeable compared to its initial state. Then, 
after some time TcoI, the wave packet has expanded a dis- 
tance 27rpniin, where Pmin is the radii of the minima of the 
sombrero. The wave packet tails will then start to over- 
lap and cause a self-interference pattern. It was shown 
in Ref. ,41] that Tcoi = 2Trv, a result which applies to the 
present model in the absence of atom-atom interaction. 
As the structure of the interference depends on the wave 
packet phases, both the dynamical phase and the Berry 
phase will be manifested in its shape. 

In order to describe the effect of the Berry phase 7, 
the dynamics rendered by the Gross-Pitaevskii equation 
pip will be compared with the dynamics of an adiabatic 
Hamiltonian possessing the same lower APS but lacking 
any Berry phase. As mentioned above, for the particular 
choice of parameters, xq = yo = Pyo = 0, Pxo = v, Ai = 1 
and ai = —02 = l/v^, the wave packet evolution takes 
place essentially on the lower APS. We therefore make 
a single surface approximation [4^ . writing an effective 
Hamiltonian as 



X y 



ad' 



(17) 



which follows from the adiabatic transformation defined 
in Eq. ([9]) and by ignoring the Berry vector potential 
^{p) — {~PvjPx)/{'2p'^) as well as the Born-Huang scalar 
potential l/(8p^) [H, 111]. Ignoring A imphes that no 
Berry phase arises when the wave packet encircles the 
CI at the origin in momentum space [4l| . On the other 
hand, the dynamical phases acquired in the two models 
are almost identical, since the effect of the Born-Huang 
potential is small for the type of evolution that is con- 
sidered here. In order to identify the consequences of the 
Berry phase on the characteristics of the two systems, it 
is therefore convenient to analyze the wave packet evolu- 
tion in the two respective models. Figure 3] depicts, for 
g = (a)-(f) and 5 7^ (g)-(l), the momentum distribu- 
tion 

MPx,Py,rf)\^ = \MPx,Py,Tf)\^ + \MPx,Py,Tf)\'^ (18) 

from our numerical propagation of the wave packet af- 
ter a time Tf — Tcoi when self-interference has been 
clearly manifested. Similar figures were presented in 
Refs. [13, Sll for different systems. The upper left and 
lower plots are obtained from the full system evolution 
governed by Eq. pT|) . while the upper right plots are 
the results achieved using the single surface approxi- 
mation P7)) . The effect of the Berry phase is evident; 



6 



With geometric phase 




Without geometric phase 







o 


IhIi 


o 






Jkm 



(b) 



« C O 



(e) <Q> Q> 






O) 











(c) 







With geometrical piiase g=-l/2 Witli geometric plrase g=l/2 




FIG. 4: Snapshots at times r = rcoi/2, 3tcoi/4, Tcoi of the 
momentum distribution Eq. (|18p after the initial wave packet 
has spread out around the minima of the sombrero poten- 
tial. Upper left plot presents the result with a Berry phase 
present and upper right one without it. Here, the dimension- 
less parameters are « = 8, t; = 0, and Tcoi = 167r. The lower 
plots (g)-(l) display the results from having a non-zero scat- 
tering amplitude 7^ 0), essentially affecting the rate of wave 
packet broadening. Again, w = 8 and Tcoi = lOvr. 



it is especially found that a node/antinode appears at 
Pa; < 0, Pj, = in the cases with/without Berry phase. 
The location of this node/antinode persists throughout 
the evolution. The total number of nodes/antinodes, on 
the other hand, changes with time, which was demon- 
strated in Ref. [41|. 

The wave packet spreading, as it evolves along the min- 
ima of the sombrero, takes place relatively far from the 
CI for the parameters of the above examples. In other 
words, the dynamics is predominantly adiabatic. How- 
ever, the concept of adiabaticity in non-linear systems is 
highly non-trivial since the superposition principle and 



orthonormality of quantum states are no longer appli- 
cable As the strength of the non-linearity is in- 
creased, the constraints for adiabatic evolution are in 
general strengthened. This has been studied in numerous 
papers, mainly focusing on various curve-crossing mod- 
els [H, 113], but also on the Berry phase [48|. In the 
current system and for the parameters considered, non- 
linearity manifests itself in the characteristic time scales 
rather than on the Berry phase effect. This is illustrated 
in Fig. m (g)-(l) where the Berry phase induced node at 
Pa; < 0, = is invariant. The number of nodes, how- 
ever, depends on the non-linearity coefficient and con- 
sequently the characteristic times are (^-dependent. This 
result may be taken as further evidence of the conjec- 
tured robustness of Berry phase effects to various kinds 
of errors and therefore of the potential utility of such 
phases in the implementation of robust quantum com- 
puters The persistence of the Berry phase effect 
due to non-linearity was also demonstrated in the model 
of Ref. [13] . We should point out that in other situations 
than the one considered in this work, namely when the 
wave packet dynamics takes place at the CI, the effect of 
the non-linearity is supposed to become important. It is 
known that non-linearity of the periodic Gross-Pitaevskii 
equation can bring about loops at the Brillouin center 
and at its edges [43|, and similar loops might be formed 
at the CI. This, however, goes beyond the scope of this 
work. 



To experimentally extract the impact of the Berry 
phase, a time-of-flight measurement can be implemented 
by switching off the trapping potential and the external 
lasers to let the condensate ballistically evolve before de- 
tecting it. In Fig. Owe show an example of such a ballistic 
expansion of the wave packet distributions \^{x,y,T)\'^ 
after it has been released from the trap. The param- 
eters are the same as in Fig. |3] (a)-(f) and the release 
time Tr = Tcoi — ISvr corresponds to Fig. 3] (c) and 
(f). The form of the momentum distributions, as shown 
in Fig. m (c) and (f), are clearly embodied in the posi- 
tion distribution. Assuming a typical trapping potential 
frequency lo /2tt ~ 40 Hz, the characteristic time scale 
is 4 ms. Hence, the time-of-fiight of Fig. [5] is around 
10 ms. This is of the same order as in common time- 
of-flight measurements of condensates in optical lattices 
[50j|. Moreover, the interaction times of Fig. 2] are of the 
order of 50-100 ms for this particular trapping potential 
frequency. Quantum phase diffusion of the condensate 
jsH , arising from quantum fluctuations beyond the mean 
field description, will most likely decrease the visibility of 
the interference pattern. However, phase coherence of 50 
ms has been observed experimentally [s^ . and at these 
time scales the Berry phase has already been manifested 
in the interference pattern. 
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FIG. 5: The position distributions after ballistic expansion. 
As in Fig. |4l the left and right plots represent the results with 
and without a Berry phase. The dimensionless parameters are 
the same as in Fig. |4] (a)-(f), and the time-of-flight duration 

is Ttof = TT. 



B. Impact of the Berry phase over long times 

The time scale Tcoi for self-interference of the wave 
packet to be established, characterizes a collapse period. 
For longer time periods, the number of wave packet max- 
ima/minima of Fig. [4] decrease to finally a single maxi- 
mum is recovered signaling a full or half revival [4l|, [5^ . 
The formation of the wave packet maxima/minima is in- 
deed highly dependent on the Berry phase. To verify 
this, we define the number of phonon excitations in the 
two modes 



_ {pi) , 



x,y. 



(19) 



This quantity is depicted in Fig. [51 where the upper three 
plots are found using the full Hamiltonian Eq. pT|) and 
the corresponding single surface results are presented in 
the lower three plots. In (a) and (d) g — 0.25, in (b) 
and (e) g = and finally in (c) and (f) g = —0.25. In- 
terestingly, at around r w 200 the population is swapped 
between the two phonon modes in the case of Berry phase 
(a)-(c) but not when the Berry phase is absent (d)-(f). At 
times r « 400 does a half revival occur [ilj. The figure 
also gives a measure on how the non-linear interaction 
affects the time scales: positive g increases while nega- 
tive g decreases the time periods over which the system 
characteristics are established. 



C. Non-Abelian manifestation 

We note that the Hamiltonian of Eq. pip can be writ- 
ten, for g = 0, as 



HcS — 



(20) 



where p = {px,Py) and A = {Ax.,Ay) = v{(Tx,(Ty). The 
gauge vector potential A has a matrix form due to the 
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FIG. 6: The average number of phonon excitations rii i = x 
(black curve) and i = y (gray curve) as a function of time 
with (a)-(c) and without (d)-(f) Berry phase. Here g = 0.25 
(a) and (d), p = (b) and (e) and g = -0.25 (c) and (f). The 
rest of the dimensionless parameters are as in Fig. |4l 



internal two-level structure of our system. The fact that 
[A^:,^^] 7^ implies that the gauge field has a non- 
Abelian character. As pointed out in Ref. [ij, encir- 
cling the CI clockwise or anti-clockwise is supposed to 
bring about different system evolution despite the po- 
lar symmetry of the system. More precisely, after one 
round-trip of the CI, the populations Pi(r) of the two 
degenerate dark states \ui) {i = 1,2) depend on the di- 
rection the wave packet traversed the minima of the som- 
brero. We demonstrate this in Fig. [7] by displaying the 
time evolution of the populations for approximately one 
round-trip clockwise (a) or anti-clockwise (b). We find 
that the progression of Pi (r) and P2 (r) is interchanged 
between clockwise or anti-clockwise wave packet evolu- 
tion, and therefore revealing the underlying non-Abelian 
structure. From the right plots, zooming in on the short 
time evolution, it is clear that the initial population is 
equally balanced between the two internal states, i.e., 
|ai|2 = |a2|2 = l/2. 
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FIG. 7: Populations Pi{T) of the dark states \ui) for clockwise 
and anti-clockwise evolution of the wave packet around the CI. 
Solid lines show the population of the state \u2} and dashed 
hnes the population of The figures to the right display 

the short time behavior, which evidence that initially both 
internal states are equally populated Pi(0) = -P2(0) — 1/2. 
The initial state are as in the previous figures {pxo = v, Ai — 1 
and ai = —02 = l/\/2) and yo = 3 (a) or yo = —3 (b). The 
other dimensionless parameters are g = and v = 4. 



is mirrored in the momentum space wave packet inter- 
ference as a Berry phase. This Berry phase has been 
computed and its effect has been identified by comparing 
the dynamics of the full system with one possessing the 
same adiabatic potential energy surface but lacking the 
Berry phase. A time-of-flight measurement would there- 
fore directly bring out the presence of such a phase. Over 
longer time periods, a periodic Berry phase dependent in- 
terchange of excitations between the x and y directions 
was demonstrated. We have also shown how the non- 
Abelian aspect of the model system is manifested in the 
wave packet dynamics in position space by propagating 
in different directions around the CI. 

It is worth mentioning that the results are not re- 
stricted to a spin-orbit coupled BEC, but can equally 
well be applied to single trapped ions [13] ■ Using the 
same laser and atomic configurations, the dynamics is 
governed by Eq. (|lip letting g ^ 0. Berry phases have 
been discussed in such systems |55i], but those settings 
are different from the one considered here. 



IV. CONCLUSIONS 

We have analyzed a spin-orbit coupled BEC system, 
emphasizing on the Berry phase and non-Abelian gauge 
structures and their importance for the dynamics. We 
have demonstrated that the non-Abelian structure as- 
sociated with the center of mass position of the atoms 



Acknowledgments 

JL acknowledge support from the MEC program 
(FIS2005-04627) and wishes to thank Prof. Stig Sten- 
holm for fruitful discussions. ES acknowledge support 
from the Swedish Research Council. 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 

Wieman, and E. A. Cornell, Science 269, 198 (1995); K. 

B. Davis, M.-O. Mewes, M. R. Andrews, N. J. vanDruten, 

D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. 

Lett. 75, 3969 (1995). 
[2] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. 

Sen, and U. Sen, Adv. Phys. 56, 243 (2007). 
[3] Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89, 130402 

(2002); T. Bourdel, L. Khaykovick, J. CubizoUes, J. 

Zhang, F. Chevy, M. Teichmann, L. Tarruell, S. J. J. 

M. F. Kokkelmans, and C. Salomon, Phys. Rev. Lett. 

93, 050401 (2004). 
[4] C. K. Law, H. Pu, and N. P. Bigelow, Phys. Rev. Lett. 

81, 5257 (1998); D. M. Stampern-Kurn, M. R. Andrews, 

A. P. Chikkatur, S. Inouye, H.-J. Meisner, J. Stenger, 

and W. Ketterle, Phys. Rev. Lett. 80, 2027 (1998). 
[5] R. Dum and M. Olshanii, Phys. Rev. Lett. 76, 1788 

(1996); P. M. Visser and G. Nienhuis, Phys. Rev. A 

57, 4581 (1998); G. Juzeliunas and P. Ohberg, Phys. 

rev. Lett. 93, 033602 (2004); P. Zhang, Y. Li, and C. 

P. Sun, Eur. Phys. J. D 36, 229 (2005); G. Juzeliunas, 

P. Ohberg, J. Ruseckas, and A. Klein, Phys. Rev. A 71, 

053614 (2005); X.-J. Liu, X. Liu, L. C. Kwek, and O. H. 

Oh, Phys. Rev. Lett. 98, 026602 (2007); M. Cheneau, S. 

P. Rath, T. Yefsah, K. J. Giinter, G. Juzeliunas, and J. 

Dalibard, Eur. Phys. Lett. 83, 60001 (2008). 



[6] K. J. Giinter, M. Cheneau, T. Yefsah, S. P. Rath, and J. 

Dahbard, Phys. Rev. A 79, 011604(R) (2009). 
[7] S. K. Dutta, B. K. Teo, and G. Raithel, Phys. Rev. Lett. 

83, 1934 (1999); Y.-J. Lin, R. L. Compton, A. R. Perry, 

W. D. Phillips, J. V. Porto, and I. B. Spielman, Phys. 

Rev. Lett. 102, 130401 (2009). 
[8] J. Ruseckas, G. Juzeliunas, P. Ohberg, and M. Fleis- 

chhauer, Phys. Rev. Lett. 95, 010404 (2005). 
[9] R. G. Unanyan, B. W. Shore, and K. Bergmann, Phys. 

Rev. A 59, 2910 (1998). 
[10] S.-L. Zhu, H. Fu, C.-J. Wu, S.-C. Zhang, and L.-M. Duan, 

Phys. Rev. Lett. 97, 240401 (2006). 
[11] A. Jacob, P. Ohberg, G. Juzeliunas, and L. Santos, Appl. 

Phys. B 89, 439 (2007). 
[12] G. Juzeliunas, J. Ruseckas, M. Lindberg, L. Santos, and 

P. Ohberg, Phys. Rev. A 77, 011803(R) (2008); M. 

Merkl, E. F. Zimmer, G. Juzeliunas, and P. Ohberg, 

Europhys. Lett. 83, 54002 (2008); G. Juzeliunas, J. 

Ruseckas, A. Jacob, L. Santos, and P. Ohberg, Phys. 

Rev. Lett. 100, 200405 (2008). 
[13] J. Y. Vaishnav and C. W. Clark, Phys. Rev. Lett. 100, 

153002 (2008). 

[14] T. D. Stanescu, C. Zhang, and V. Galitski, Phys. Rev. 

Lett. 99, 110403 (2007). 
[15] J. Y. Vaishnav, J. Ruseckas, C. W. Clark, and G. Juzelu- 

nas, Phys. Rev. Lett. 101, 265302 (2008). 



9 



[16] L.-H. Lu and Y.-Q. Li, Phys. Rev. A 76, 023410 (2007). 
[17] T. D. Stanescu, B. Anderson, and V. Galitski, Phys. Rev. 

A 78, 023616 (2008). 
[18] A. Jacob, P. Ohberg, G. Juzeliunas, and L. Santos, New 

J. Phys. 10, 045022 (2008). 
[19] V. Pietila and M. Mottonen, Phys. Rev. Lett. 102, 

080403 (2009). 

[20] M. Baer, Beyond Born-Oppenheimer (Wiley, New York, 
2006). 

[21] C. A. Mead, Rev. Mod. Phys. 64, 51 (1992). 

[22] C. A. Mead and D. G. Truhlar, J. Chem. Phys. 70, 2284 

(1979); C. A. Mead, Chem. Phys. 49, 23 (1980). 
[23] M. V. Berry, Proc. R. Soc. London Ser. A 392, 45 (1984). 
[24] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005) ; Y. B. Zhang, Y. W. Tan, H. L. Stormer, and 
P. Kim, Nature 438, 201 (2005); K. S. Novoselov, A. 
K. Geim, S. V. Morozon, D. Jiang, M. L Katsnelson, L 
V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 
438, 197 (2005). 

[25] A. D. Dudarev, R. B. Diener, L Carusotto, and Q. Niu, 
Phys. Rev. Lett. 92, 153005 (2004); S. L. Zhu, B. G. 
Wang, and L. M. Duan, Phys. Rev. Lett. 98, 260402 
(2007). 

[26] H. C. Longuet-Higgins, U. Opik, M. H. L. Pryce, and R. 

A. Sack, Proc. R. Soc. London Ser. A 244, 1 (1958). 
[27] D. R. Yarkony, Rev. Mod. Phys. 68, 985 (1996). 
[28] A. Bohm, A. Mustafazadeh, H. Kuizumi, Q. Niu, and 

J. W. Zwanziger, Geometric Phase in Quantum Systems 

(Springer, Berlin, 2003). 
[29] L Fuentes-Guridi, J. Pachos, S. Bose, V. Vedral, and S. 

Choi, Phys. Rev. A 66, 022102 (2002); E. I. Duzzioni, L. 

Sanz, S. S. Mizrahi, and M. H. Y. Moussa, Phys. Rev. A 

75, 032113 (2007). 
[30] Z. D. Chen, J. Q. Liang, S. Q. Shen, and W. F. Xie, 

Phys. Rev. A 69, 023611 (2004). 
[31] J. Samuel and R. Bhandari, Phys. Rev. Lett. 60, 2339 

(1988); N. Mukunda and R. Simon, Ann. Phys. (N.Y.) 

228, 205 (1993). 
[32] J. Larson and S. Stenholm, Phys. Rev. A 73, 033805 

(2006) . 

[33] Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 
(1984). 

[34] S.-R. E. Yang and N. Y. Hwang, Phys. Rev. B 73, 125330 
(2006); P. Pietilainen and T. Chakroborty, Phys. Rev. B 
73, 155315 (2006). 

[35] H. Tiitiinciiler, R. Koc, and E. Olgar, J. Phys. A: Math. 
Gen. 37 11431 (2004). 

[36] H. Feshbach, Ann. Phys. 5, 337 (1958); E. Tesinga, B. 
J. Verhaar, and H. T. C. Stoof, Phys. Rev. A 47, 4114 
(1993); C. Chin, R. Grimm, P. Julienne, and E. Tisinga, 
larXiv:0812. 14961 



[37] M. D. Flei, J. A. Fleck, and A. Steiger, J. Comp. Phys. 
47, 412 (1982). 

[38] D. M. Tong, E. Sjoqvist, L. C. Kwek, and C. H. Oh, 

Phys. Rev. Lett. 93, 080405 (2004). 
[39] E. Sjoqvist and M. Hedstrom, Phys. Rev. A 56, 3417 

(1997); R. Englman, A. Yahalom, and M. Baer, Eur. 

Phys. J. D 8 1 (2000). 
[40] J. Schon and H. Koppel, J. Chem. Phys. 103, 9292 

(1995) ; J. Schon and H. Koppel, J. Chem. Phys. 108, 
1503 (1998). 

[41] J. Larson, Phys. Rev. A 78, 033833 (2008). 

[42] R. Baer, D. M. Charutz, R. Kosloff, and M. Baer, J. 
Chem. Phys. 105, 9141 (1991); S. Adhikari and G. D. 
Billing, J. Chem. Phys. Ill, 40 (1999). 

[43] M. Born and K. Huang, Dynamical theory of crystal lat- 
tices (Clarendon, Oxford, 1954), pp. 406-407. 

[44] M.V. Berry and R. Lim, J. Phys. A 23, L655 (1990). 

[45] J. Liu, B. Wu, and Q. Niu, Phys. Rev. Lett. 90, 170404 
(2003). 

[46] B. Wu and Q. Niu, Phys. Rev. A 61, 023402 (2000); V. 
A. Brazhnyi and V. V. Konotop, Mod. Phys. Lett. B 18, 
627 (2004). 

[47] M. Machholm, C. J. Pethick, and H Smith, Phys. Rev. 
A 67, 053613 (2003). 

[48] B. Wu, J. Liu, and Q. Niu, Phys. Rev. Lett. 94, 140402 
(2005); X. X. Yi, X. L. Huang, and W. Wang, Phys. Rev. 
A 77, 052115 (2008). 

[49] J. A. Jones, V. Vedral, A. Ekert, and G. Castagnoli, 
Nature 403, 869 (2000); W. Xiang-Bin and M. Keiji, 
Phys. Rev. Lett. 87, 097901 (2001); D. Leibfried, B. De- 
Marco, V. Meyer, D. Lucas, M. Barrett, J. Britton, W. 
M. Itano, B. Jelenkovic, C. Langer, T. Rosenband, and 
D.J. Wineland, Nature 422, 412 (2003). 

[50] F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, 
and I. Bloch, Phys. Rev. Lett. 95, 050404 (2005). 

[51] M. Lewenstein and L. You, Phys. Rev. Lett. 77, 3489 

(1996) . 

[52] D. S. Hall, M. R. Matthews, C. E. Wieman, and E. A. 

CorneU, Phys. Rev. Lett. 81, 1543 (1998); M. Albiez, R. 

Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M. K. 

Oberthaler, Phys. Rev. Lett. 95, 010402 (2005). 
[53] R. W. Robinett, Phys. Rep. 392, 1 (2004); E. Romera 

and F. de los Santos, Phys. Rev. Lett. 99, 263601 (2007). 
[54] D. Leibfried, R. Blatt, C. Monroe, and D. J. Wineland, 

Rev. Mod. Phys. 75, 281 (2003). 
[55] 1. Fuentes-Guridi, S. Bose, and V. Vedral, Phys. Rev. 

Lett. 85, 5018 (2000); X.-D. Zhang, Z. D. Wang, L.-B. 

Hu, Z.-M. Zhang, and S.-L. Zhu, New J. Phys. 10, 043031 

(2008). 



